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ABSTRACT 

To measure the properties of both components of the RS CVn binary a Geminorum (u Gem), we 
directly detect the faint companion, measure the orbit, obtain model-independent masses and evolu¬ 
tionary histories, detect ellipsoidal variations of the primary caused by the gravity of the companion, 
and measure gravity darkening. We detect the companion with interferometric observations obtained 
with the Michigan InfraRed Combiner (MIRC) at Georgia State University’s Center for High Angular 
Resolution Astronomy (CHARA) Array with a primary-to-secondary iJ-band flux ratio of 270 ± 70. 
A radial velocity curve of the companion was obtained with spectra from the Tillinghast Reflector 
Echelle Spectrograph (TRES) on the 1.5-m Tillinghast Reflector at Fred Lawrence Whipple Observa¬ 
tory (FLWO). We additionally use new observations from the Tennessee State University Automated 
Spectroscopic and Photometric Telescopes (AST and APT, respectively). From our orbit, we deter¬ 
mine model-independent masses of the components (Mi = 1.28±0.07 Mq, M 2 = 0.73±0.03 Mq), and 
estimate a system age of 5 =F 1 Gyr. An average of the 27-year APT light curve of cr Gem folded over 
the orbital period {P = 19.6027 ± 0.0005 days) reveals a quasi-sinusoidal signature, which has previ¬ 
ously been attributed to active longitudes 180° apart on the surface of a Gem. With the component 
masses, diameters, and orbit, we find that the predicted light curve for ellipsoidal variations due to 
the primary star partially filling its Roche lobe potential matches well with the observed average light 
curve, offering a compelling alternative explanation to the active longitudes hypothesis. Measuring 
gravity darkening from the light curve gives /3 < 0.1, a value slightly lower than that expected from 
recent theory. 

Subject headings: binaries: close - stars: activity - stars: imaging - stars: individual (cr Geminorum) 
- stars: variables: general 


1. INTRODUCTION 

RS Ganum Venaticorum (RS GVn) stars are spot¬ 
ted, active binary systems exhibiting photometric and 
Ga H and K variability (jHalll Il976fl . Often tidally- 
locked, these systems are composed of an evolved pri¬ 
mary star ( giant or subgiant) and a subgiant o r dwarf 
companion i|Berdvuginal l2005t iStrassmeied I 200l). With 
active binaries, not only is there potential to deter¬ 
mine the component masses and system evolutionary 
history but also to understand the magnetic field in¬ 
teractions through active longitudes, particular longi- 
tudes 180° apart with persistent, long-lived sta rspots 
(|Berdvugina fc Tuomineniri998l : lBerdvugin^l2005[l . 

Observing the magnetic phenomena of rapidly-rotating 
evolved stars also sheds light on the magnetic activity 
of rapidly-rotating young stars, such as T Tauri stars. 
Both T Tauri and RS GVn systems have starspots anal¬ 
ogous to sunspots—cool starspots resulting from stifled 
convection in the outer layers of the stars due to strong 
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magnetic fields (|Petrovll2003l: lBerdvuginall2005ll . 

cr Geminorum (cr Gem, HD 62044, HIP 37629, 
HR 2973) is an RS GVn system known to exhibit 
starspots, often ascribed to “activ e longitudes” (e.g., 
iHall h ahl 119771 : iHenrv et al.l I1995D . The system has 
been characterized as a si ngle-lined spectroscopi c binary 
(iHerbig fc SDaldindll955t) with a KlHl primary (|Romanl 
li952ll . The orbital period of cr Gem is slightly longer 
than the primary star’s rotation period derived from the 
fastest rotating spots (Pnrh = 19.60 days, Pmt min = 
19.47 days; iKaiatLriet akfeoTl . 

Because of its large starspots, a Gem is a fre- 
quent target for unders t andiii g starspot evolution. 
lEberhard fc SchwarzschiO (jl913f l first reported cr Gem 
as active and potentially spotted due to fluctuations in 
the Ga H and K line s as the star rotated. Decades 
later, iHall et ahl (|1977f) identified photometric variations 
suggesting starspots (AU ^ 0.07). Initial models of 
the surface of cr Gem often showed the surface with 
two starsp ots oriented on o pposite sides of the pri- 
mary star ((Fried et al.l 1198311 . iBerdvugina fc Tuominenl 
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(|1998fl emphasize that, due to tidal locking, the starspots 
are located such that one spot constantly faces the 
companion and the other spot is 180° offset. The 
majority of spot models applied to light curves of 

a Ge m _ con sist of two spots on a spherica l star 

(lEked 119801: iStrassineier et al.l 119881 lOlah et all 119891 
Henry et al.i 11995: iJetsul 119961 : iPadmakar fc PandevI 
19991: iKaiatkari et al.ll20l4l . Doppler images have sug¬ 
gested the su rface is covered with a larger numbe r of 
smaller spots (jHatzeslll99,‘ll : lKdvari et al.l 1200 IL[2?)T^ . 

To understand the binary system, we present our anal¬ 
ysis of the first detections of the companion in our inter¬ 
ferometric and radial velocity data sets, as well as photo¬ 
metric data. In Section 2, we describe the observations 
for our data sets. In Section 3, we discuss our analy¬ 
sis of the data sets, including the first astrometric and 
spectroscopic detections of the companion star and or¬ 
bital parameters. In Section 4, we present evolutionary 
constraints and a Hertzsprung-Russell (H-R) diagram. 
In Section 5, we discuss our analysis of the photomet¬ 
ric data set, including detected ellipsoidal variations and 
measured gravity darkening. In Section 6, we present the 
conclusions of our study of a Gem. 

2. OBSERVATIONS 

2.1. Interferometry 

We obtained interferometric data with Georgia State 
University’s Genter for High-Angular Resolution Astron¬ 
omy (GHARA) Array. The GHARA Array is a Y-shaped 
array of six 1-m class telescopes with non-redundant 
baselines varying from 34- to 33 1-m located at Mount 
Wilso n Observatory, Galifornia (|ten Brummelaar et al.l 
120051) . Using all six tel escopes and the Michigan In- 
fraRed Gombiner ('MIR.Gi rMoimier et al.l[20(i^ l2006| ) . we 
obtained iJ-band (1.5 — 1.8 /rm) data (eight channels 
across the photometric band with A/AA ^ 40) on UT 
2011 Nov 9 and Dec 7, 8, 9; 2012 Nov 7, 8, 21, 22, 24, 25 
and Dec 4, 5. 

We made detections of the companion in the data from 
UT 2011 Dec 8; 2012 Nov 7, 8, 24, and 25. The re¬ 
maining nights of observation had insufficient uv cov¬ 
erage due to poor seeing or short observation lengths, 
leaving the companion undetected. We reduced and cal- 
ibrated these data w i th the standard MIRG p i peline (see 
Monnier et al.l 120071 120121: iZhao et HI 120091 iGhe et all 
20 111 for pipeline details). We used at least one calibra¬ 
tion star for each night of data (see Table [IJ. 

2.2. Radial Velocity 

To constrain the spectroscopic orbit for a Gem, we 
utilized three independent sets of radial-velocity data: 
two sets of single-lined velocities for the primary, and a 
new set of double-lined velocities for both components of 
the binary. 

One set of the radial veloci ty measurements for th e 
primary star was published in iMassarotti et al] (1200811 . 
These 39 data points were obtained with two identical 
GfA Digital Speedometers (jLathamI 119921 1 on the 1.5- 
m Wyeth Reflector (Oak Ridge Observatory) and 1.5-m 
Tillinghast Reflector (Fred Lawrence Whipple Observa¬ 
tory) telescopes (2003 December 30 — 2007 June 5). 

From 2012 October 1 — 2015 January 9, using the Till¬ 
inghast telescope with the Tillinghast Reflector Fchelle 


Spectrograph (TRFS; iFiireszl [200811 . we were able to 
make fifteen detections of the secondary spectra for the 
first time. Along with sixteen new primary star measur- 
ments, these new radial velocities are presented in Table 
[21 We add 0.14 kms“^to these sets of radial velocities 
to account fo r these data being re ported on the GfA na¬ 
tive system (iStefanik et al.l [19^ note the correction is 
inaccurately stated as a subtraction in this reference). 
For details of these observations and data analysis, see 
Appendix A. 

The additional radial velocity data set consists of 43 
spectrograms of the primary star of a Gem taken between 
2009 January 12 — 2014 December 1 with the Tennessee 
State University 2-m automatic spectroscopic telescope 
(AST), fiber-fed echelle spectrograph, and a GGD de- 
tector at Fairborn Observatory, Arizona (see Table |3l 
iFaton fc WilliamsonI 12004 1200711 . At first, the detector 
was a 2048 x 4096 SITe ST-002 A GGD with 15 /rm pix¬ 
els. [Eaton fc WilliamsonI (1200711 discussed the reduction 
of the raw spectra and wavelength calibration. Those 
echelle spectrograms have 21 orders that cover the wave¬ 
length range 4920-7100 A with an average resolution of 
0.17 A, corresponding to a resolving power of 35000 at 
6000 A. Those spectra have a typical signal-to-noise value 
of 30. 

In the summer of 2011 the AST SITe GGD and its de- 
war were retired and replaced with a Fairchild 486 CCD, 
a 4096 X 4096 array of 15 /im pixel, that is housed in a 
new dewar. With the new CCD the wavelength coverage 
ranged from 3800 to 8600 A. The resolution was reduced 
slightly to 0.24 A or a resolving power of 25000 at 6000 A. 
These more recent spectra have signal-to-noise ratios of 
a bout 70. 

iFekel et ahl (|2n09[l provided an extensive general de¬ 
scription of velocity measurement of the Fairborn AST 
spectra. In the case of a Gem, we measured a subset of 
63 lines from our solar-type star line list that covers the 
4920-7120 A region. Because the lines of a Gem have 
significant rotational broadening, we fit the individual 
lines with a rotational broadening function. The Fair¬ 
born velocities are on an absolute scale. A comparison of 
our unpublished measurements of several lAU sta ,ndard 
stars with those determined bv iScarfe et al.l (|1990tl indi¬ 
cates that the Fairborn Observatory velocities from the 
SITe CCD have a small zero-point offset of —0.3 kms“^. 
Velocities from the Fairchild CCD spectra have a slightly 
larger zero-point off set of —0.6 kms“^relative to those of 
iScarfe et'HI (|1990(1 . Thus, in Table |3| we corrected our 
measured velocities by either 0.3 or 0.6 kms“^, depend¬ 
ing on which detector was used. 

2.3. Photometry 

We used differential photometry of cr Gem and a com¬ 
parison star from the Tennessee State University T3 0.4- 
m Automated Photometric Telescope (APT) located at 
Fairborn Observatory, Arizona. For deta ils on the ob ¬ 
serv ational procedure and photometers see lHenr^ (Il999f l 
and iFekel et al.l ()2005l l. 

The differential Johnson B and V light curves cover 
1987 November 21 — 2015 March 13 (see Table HI and 
Figure H]). Subse ts o f these data we r e ana lyzed by 
iHenrv et al.l (I1995D and IKaiatkari et HI (1201411 . For the 
first time, we make the full set of T3 APT photometry 
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TABLE 1 

Calibrators for a Geminorum 


Calibrator Name 

Calibrator Size (mas) 

Source 


UT Date of Observation 

HD 37329 

HD 50019 [B Gem) 
HD 63138 

HD 69897 (x Gnc) 

0.71 ±0.05 

0.81 ±0.06 

0.65 ±0.04 

0.73 ± 0.05 

Bonneau et al. ('2006') 
Bonneau et al. ^556') 
MiHC calibration 
Bonneau et al. f2006') 

2012 Nov 8 

2012 Nov 7, 8, 25 

2011 Dec 8; 2012 Nov 8 
2012 Nov 7, 24, 25 


TABLE 2 

Radial Velocity Data of a Gem (CfA) 


HJD -2400000 

Primary (kms 

Secondary (km s ^) 

56202.0199 

10.66 

103.44 

56230.0436 

77.98 

-9.80 

57002.9567 

14.17 

94.68 

57003.9301 

9.71 

101.50 

57014.8600 

79.02 

-12.11 

57015.8427 

75.96 

-10.77 

57018.8519 

51.08 

28.91 

57019.9181 

38.76 

60.57 

57020.9452 

27.84 

78.02 

57021.8962 

19.05 

84.20 

57024.9812 

9.10 

104.10 

57025.9317 

12.59 

99.41 

57026.8891 

18.88 

85.69 

57028.9057 

38.58 


57029.8670 

49.81 

27.28 

57031.8806 

69.76 

-2.88 


Note. — Errors on the primary radial velocities are 0.84 
kms“^. Errors on the secondary radial velocities are 3.8 kms“^. 
These were then scaled for our orbit fit to have a total = 1.00. 
Note these radial velocities are on the n ative CfA system. W e 
added 0.14 kms~^for use in our analysis iStefanik et al.1l 19990 . 



Fig. 1. — Johnson B and V differential magnitudes of g Gem 
acquired over 28 observing seasons from 1987 — 2015 with the T3 
0.4-meter APT at Fairborn Observatory in southern Arizona. 

of a Gem available in electronic format in Table H] 

3. ORBITAL ELEMENTS 

In order to derive the astrometric orbit of a Gem, we 
searched for the companion with model fitting. We mod¬ 
eled the system with the resolved primary star and an 
unresolved secondary. We allowed the primary radius 
along the major axis, primary major-to-minor axis ratio, 
primary major axis position angle, primary-to-secondary 
flux ratio, and secondary position to vary. During the fit¬ 


TABLE 3 

Radial Velocity Data of a Gem 
(AST/TSU) 


HJD -2400000 

Primary (kms 

54843.9716 

38.0 

54844.8597 

28.4 

54845.6849 

20.5 

54846.6455 

14.0 

54847.6848 

9.2 

54848.7712 

9.2 

54849.7451 

12.1 

54850.6645 

17.6 

54856.9772 

75.7 

54859.7922 

73.1 

54860.6910 

66.7 

54861.6447 

58.7 

54862.6829 

48.3 

54863.6434 

37.3 

54864.6443 

26.7 

54865.6596 

17.8 

54866.6401 

11.7 

54867.6456 

9.1 

54868.9572 

10.3 

54869.6594 

13.8 

54870.6346 

21.1 

54975.6558 

77.9 

54976.6372 

76.4 

54981.6497 

33.3 

54982.6541 

23.0 

54983.6544 

15.2 

54984.6546 

9.9 

54985.6548 

8.5 

54986.6547 

10.7 

55060.9627 

23.5 

55061.9642 

16.0 

55066.9516 

23.3 

55069.9622 

54.3 

55070.9596 

64.0 

55071.9572 

71.9 

55072.9375 

76.6 

56984.7820 

9.4 

56987.7638 

19.7 

56988.7240 

28.0 

56989.7317 

38.9 

56990.7354 

50.0 

56991.7355 

60.9 

56992.8661 

70.8 


Note. — Errors on the primary ra¬ 
dial velocities are 0.3 kms”^. These data 
were then scaled for our orbit fit to have 
a total x^ — 1-00. 

ting, we weighted the data such that the separate observ¬ 
ables (squared visibilities, closure phases, and triple am¬ 
plitudes) contributed to the final with equal weight. 
The parameter errors for the primary star size and the 
primary-to-secondary flux ratio were based on the epoch- 
to-epoch variation, while the relative positional error of 
the secondary compared to the primary were based on 
the residuals to the orbit fit (see discussion on orbit fit¬ 
ting). 

The coordinates of the detections on five nights (UT 
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TABLE 4 

Johnson B and V Differential Magnitudes of a Gem from the T3 APT 


HJD -2400000 

AB (cr Gem - HD 60318“) 

AV (a Gem - HD 60318) 

AB {v Gem*" - HD 60318) 

AH {v Gem - HD 60318) 

47121.0481 

-1.006 

-1.123 



47122.0465 

-1.010 

-1.135 



47125.0364 

-1.030 

-1.125 

-0.766 

-1.303 

47126.0364 

-1.011 

-1.128 

-0.757 


47128.0310 

-0.975 


-0.759 



Note. — Table^is published in its entirety in the electronic edition. A portion is shown here for guidance regarding its form and content. 
“ HD 60318, V = 5.33, B - V = 1.01 
^ V Gem (HD 60522, 69 Gem), V — 4.06, B — V — 1.54 


2011 Dec 8; 2012 Nov 7, 8, 24, and 25 are listed in Ta¬ 
ble [S]). The i7-band flux ratio for the primary star to 
the secondary is 270 ± 70. In addition to detecting the 
secondary star, we measured the uniform disk diameter 
of the primary to be 0 ud,i = 2.335 ± 0.007 mas (limb- 
darkened disk diameter 6 *ld,i = 2.417 ± 0.007 mas) with 
a major-to-minor axis ratio of 1.02 ± 0.03. Our measure¬ 
ments are slightly larger than those in the CHARM2 cat¬ 
alog (uniform disk diameter of duD,i = 2.18 ± 0.05 mas, 
limb-darkened disk diameter of ^ld,! = 2.31 ± 0.05 mas; 
iRichichi et al.l[2?)?)^ . 

To determine the binary orbit, we simultaneously fit 
our interferometric and radial velocity data with Monte 
Carlo realizations. The five interferometric points are as 
described above, and we present the scaled error bars of 
the major and minor axis in Table [5] to give our fit a 
total = 1.00. For the radial velocity data we com¬ 
bine the iMassarotti etldl (j2008l adding 0.14 kms“^to 
account for the values reported on the CfA native sys¬ 
tem) , new CfA data, and the AST data to fit simultane¬ 
ously with the astrometry. The radial velocity errors are 
similarly scaled (rmscfA,i = 0.84 kms“^, rmsAST,i = 0.3 
kms“^, rmscfA,2 = 3.8 kms“^). 

Using the complete radial velocity data sets, we 
find an eccentricity of e = 0.014 ± 0.004, consis¬ 
tent w it h slightly eccen t ric or bits rep orted bv iHaroei 
( 1935D . iPourbaix et al.l (200411. and IMassarotti et al. 
( 200^. However, Luyten ( 1936f l. iBatten et al.l ( ~78ll . 
and iDiimmler et al.l (jl997^ reported a circular orbit. 
To investigate this discrepancy, we used the APT light 
curve to eliminate the primary star’s radial velocity data 
that were obtained when a Gem presented starspots 
(AU > O. C)4), as these could cau se shifts in the veloc¬ 
ities fe.g.. iSaar fc Donahuelll99'^ . The remaining pri¬ 
mary star radial velocity data obtain ed when a Gem 
did no t exhibit large starspots from the IMassarotti et al.l 
(j2008fl /CfA data set span 2006 December 6 — 2007 June 
5, and those from the AST data set span 2009 January 
12 — Jun 4. Us ing the primary s t ar’s tr uncated data set 
with 42% of the lMassarotti et '5l (|2008ll /CfA and 33% of 
the AST epochs removed, we find the orbit is consistent 
with a circular orbit, e = 0.002 ± 0.002, and we adopt a 
circular orbit for the rest of this paper. 

Requiring eccentricity e = 0 and the argument of pe- 
riastron for the primary w = 0° the simultaneous Monte 
Carlo realizations gave the orbital parameters and their 
1 — a errors listed in Table The visual orbit is il¬ 
lustrated in Figure [21 and the radial velocity curve is 
pr esented i n Figu re |3l We use the conventions presented 
bv iHeinti (119781 1. where the argument of periastron, w, 
and the time of nodal passage (maximum recessional ve¬ 


locity), To, are defined by the primary star’s orbit. The 
ascending node, fl, is independent of definition, being 
equivalent with respect to either the primary or sec¬ 
ondary star. 

Our orbital parallax, tt = 25.8 ± 0.4 mas can be 
comp ared with t he Hipparcos parallax of 26.68 ± 0.79 
mas (If^ [1991). As an unresolved binary with a 
variable component, a Gem does not exhibit the pho¬ 
tocenter shifts found to be troublesome f or measur¬ 
ing binary system parallax w ith Hipparcos (lESAl[l9^ 
iHalbwachs fc Poiirbaixl 1200511 . Assuming that the sec¬ 
ondary is negligibly bright, the semi-major axis of the 
photocentric orbit of the primary is at most 1.71 mas 
wide, which is at the limit of detectability ([Pourbaixl 
I 2 OO 2 II for Hipparcos. Gombining Hipparcos data and our 
visual orbit, the parallax is 26.4±0.8 mas, consistent with 
our orbital parallax. For our subsequent analysis, we 
adopt our higher-precision orbital parallax, tt = 25.8±0.4 
mas. 

With a circular orbit and JUrb Trot (e.g., 

iKaiatkari et al.ll2014ll . we expect a Gem to have aligned 
rotational and orbital axes. Given our orbital and stel¬ 
lar parameters, we can calculate the obliquity of the 
system. Comparing our calculated value of v sin i = 
(27ri?i/Porb) X usinf = 24.8 ± 0.4 kms“^with the ob¬ 
servational rotational velocity of vsini = 26.7 ± 0.5 
kms“^(from the TRES spectra), we find that the cal¬ 
culation is smaller than the observational value. This 
discrepancy could be attributed to the estimate of mir- 
croturbulence or the presence of the large spot structures 
on the surface of a Gem during the TRES observations 
instead of a small, non-zero obliquity. 

4. MASSES AND HERTZSPRUNG-RUSSELL DIAGRAM 

Using our complete orbital fit, we obtain model- 
independent masses Mi = 1.28 ± 0.07 Mq and M 2 = 
0.73 ± 0.03 Mq. With the stellar parameters of the pri¬ 
mary star (including Teffp = 4530 ± 60 K, see Table [6]) 
and the primary-to-secondary JJ-band flux ratio detected 
using the CHARA/MIRC data (270 ±70), we are able to 
constrain the parameters (luminosity, temperature, and 
radius) of the secondary star. We use the flux ratio a nd 
NextGen stellar atmospheres ([Hauschildt et al.lll999ll to 
constrain the stellar flux to calculate a range of luminosi¬ 
ties for reasonable effective temperatures (4000 — 4700 K) 
for a 0.73 ± 0.03 Mq main sequence star (see Figured]). 
We obtain a range of luminosities (0.11 — 0.15 Lq) and 
radii (0.70 — 0.59 Rq). We note that our analysis predicts 
a primary-to-secondary Johnson U-band flux ratio of 290 
assuming Teff _2 = 4500 K), which is not in agreement 
with the flux ratio given by the spectroscopic 519 nm 
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TABLE 5 

Detections for the companion of o- Geminorum with respect to the primary 


UT Date 

JD -2400000 

Separation 

(mas) 

Position 
Angle (°)“ 

Error Ellipse 
Major Axis (mas)^ 

Error Ellipse 
Minor Axis (mas)^ 

Error Ellipse 
Position Angle (°)“ 

Reduced 

2011 December 08 

55903.95 

2.83 

19.1 

0.30 

0.09 

80 

4.4 

2012 November 07 

56238.97 

4.32 

8.6 

0.04 

0.03 

280 

2.8 

2012 November 08 

56239.86 

4.68 

359.7 

0.13 

0.06 

300 

2.1 

2012 November 24 

56256.00 

2.03 

39.0 

0.08 

0.06 

30 

1.5 

2012 November 25 

56256.95 

3.12 

21.6 

0.05 

0.04 

320 

1.7 


Note. — These detections give an H^-band (1.5 — 1.8 flux ratio for cr Gem primary to secondary of 270 ± 70. The uniform disk fit for the 
primary star is 0uD,i = 2.335 i 0.006 mas (limb-darkened disk diameter 0 ld,i = 2.417 ± 0.006 mas) with a 1.02 ± 0.03 major-to-minor axis ratio. 
“East of North 

^Scaled error bars to ensure a total = 1.00 as described in Section 3.2. 


TABLE 6 

Orbital and Stellar Parameters of <j Gem 


Measured Parameters 

Value 

semi-major axis, a (mas) 

eccentricity, e 

inclination, i {°) 

argument of periastron, u (°)“ 

ascending node, H (°) 

period, Poj-b (days) 

time of nodal passage. To (HJD)^ 

velocity semi-amplitude, Ki (kms“^) 

velocity semi-amplitude, K 2 (kms“^) 

system velocity, 7 (kms“^) 

uniform disk diameter, ^uD,i (mas) 

limb-darkened disk diameter, ^ld,i (mas)*^ 

primary major-to-minor axis ratio 

P-band flux ratio, primary to secondary 
orbital parallax, tt (mas) 
distance, d (pc) 

4.63 ± 0.04 

0 

107.7 ±0.8 

0 

1.2 ± 0.8 
19.6027 ± 0.0005 
2453583.98 ± 0.03 
34.62 ±0.08 

60 ±2 

43.41 ± 0.08 
2.335 ±0.007 
2.417 ±0.007 

1.02 ±0.03 

270 ± 70 

25.8 ±0.4 

38.8 ±0.6 

Derived Parameters 

average primary radius, Ri (Rq)^ 
primary luminosity, Li (Tq) 
primary surface gravity, \oggi (cm/s^) 
primary mass, Mi (Mq) 
secondary mass, M 2 (Mq) 
system age (Gyr) 

10.1 ±0.4 

39 ±2 

2.54 ±0.02 

1.28 ±0.07 

0.73 ±0.03 

5±1 

Literature Parameters 

primary effective temperature, (R)^ 

primary metallicity (iron), Fe/H^ 

4530 ± 60 

0.0 

Note. — “Radial velocity convention for primary with respect to 
the center of mass. 

^Time of maximum recessional velocity of the primary star. 

“We applied a 3.5% correction from uniform to limb-darkened disk 
diameter. This is equivalent to a limb-darkening coefficient q = 0.27. 


‘^Using limb-darkened disk diameter. 

^Teff'i is an av erage of te mperatures given _ by 

Glgbock^^^tawikow^^ lll9 79p: _ [Poe^^Eat^ _ 

Stawilcowsk^^^lgbo^ (1^94 b: lo Neal et aTT^^g^ : iKovari et al.| 
j200TIT 7 lMassarott^t^U 1200^ . The 1 — cr error is the standard 



Fig. 2.— Visual orbit for the prototypical RS CVn system 
a Gem with our observed stellar primary radius (thick black line, 
cr Gem A) and our dates of companion detection and their loca¬ 
tions on the orbit (black error ellipses). The predicted radius of 
the companion star, cr Gem B, is plotted for scale with the small 
thick black circle. The orbits of fifty Monte Carlo realizations are 
presented as the light gray orbits. Black lines connect the center 
of the detection error ellipse to the expected point in the best-fit 
orbit, which is overlaid in black (given in Table with 1 — cr er¬ 
rors). At the southernmost point in the orbit, the secondary star is 
moving toward the observer. Note: axis units are milliarcseconds 
(mas) with north upwards and east to the left. 


deviation of these values. _ 

■^[Fe/H] — —0.02 ||MaIlik||1998f) ; approximated as [Fe/H]— 0.00. 

light ratio 70 primary-to-secondary; see Appendix A 
for details on this measurement). In order for our flux 
ratios to be in agreement, the secondary star would have 
7eff,2 = 6400 K, which is not consistent with the spectro¬ 
scopic observations, nor with a main-sequence star given 
the location on the H-R diagram. We cannot rule out the 
effect of starspots on the discrepant flux ratios as these 
were not accounted for when interferometrically detect¬ 
ing the companion and the spot features present dur¬ 
ing the interferometric and spectroscopic observations 
diffe r as evidenc e d in the APT light curve. Add ition- 
ally, iPrato et'afl (1200211 and iLehmAn et al.l (1201311 also 
reported discrepancies between TODCOR-reported flux 


ratios and their expected values. Therefore, we use only 
the i7-band flux ratio. 

We plot the location of the components of a Gem 
on an H-R diagram, as well as the corresponding evo¬ 
lutionary tracks. We use Dartmouth stellar evolution 
tracks ('Fe/H = 0.(3 g/F e = 0.0, PHOENIX-based 
models; iDotter et al.l l200§l for the interpolated model 
masses (Mi,model = 1.28 ±0.07 M©, M 2 ,model = 0.73 ± 
0.03 Mq). Our primary falls nearly on the 1.28 Mq 
evolutionar y track with an estimated temperature of 


4530±60K (IGlebocki & Stawikowsk 

1979HPoe & EatonI 

19851: IStawikowski & Glebockil 1991 

O’Neal et al.lll996l 

Kovari et al.ll200lHMassarotti et al. 

200811. The ransre of 

locations for the secondary on the H-R diagram passes 
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Fig. 3. — Radial velocity curves of a Gem. The filled di- 
amonds represent our sample of measured observations from 
IMassarotti et aP 1)20081) /CfA. and the filled circles are the AST 
observations. Both data sets are restricted to those data points 
obtained with no starspots present (see Section 3). 1 — a errors in 
velocity are presented unless the error is smaller than the diamonds 
and circles. The radial velocity curves of fifty Monte Carlo realiza¬ 
tions are presented as the light gray orbits. The radial velocity for 
the best orbital parameters is overlaid in black. Similarly, the open 
diamonds represent CfA radial velocity data for the secondary star 
with 1 — cr error bars. The light gray orbits are fifty Monte Carlo 
realizations with the best orbital parameters overlaid in black. See 
Table [ 6 ] for orbital parameters with 1 — a errors. 



Teff{K) 


Fig. 4.— H-R diagram for cr Gem. The dashed and dot-dashed 
lines are the main sequence and post-main sequence evolutionary 
tracks for 1.28 M^t) a nd 0.73 Mq stars with [Fe/H] 0.0, respec¬ 
tively JDotter et al.|[^08l) . The gray regions represent our 1 — cr 
mass errors (Mi = 1.28 zb 0.07 Mq;M 2 = 0.73 ib 0.03 Mq) with 
the solid black line representing the zero age main sequence. Th e 
dotted line is a 5 Gyr isochrone (PHOENIX: [Potter et ^120081) . 
The measured location of the primary with 1 — cr errors is indi¬ 
cated by the plus sign. The region where the companion could 
be located given our flux ratio and reasonable temperature esti¬ 
mates is indicated with the long-dashed line (with 1 — cr errors in 
luminosity). 


through the main sequence for a star of 0.73 Mq. We 
find an age of the system of 5 =F 1 Gyr. Based upon the 
masses and age of the stars, we suggest that the primary 
star is an evolved late F-type star that is now a K giant. 
The secondary star is a main-sequence early K star. 

5. ELLIPSOIDAL VARIATIONS AND GRAVITY 
DARKENING 

iHenrv et al.l (119951) and iKaiatkari et al.l (I2014D previ¬ 
ously published subsets of the APT light curve data for 
starspot modeling and measuring differential rotation. 
Both studies emphasized the presence of active longi¬ 
tudes on opposite sides of a Gem to explain the quasi- 
sinusoidal variation appearing at half of the orbital pe¬ 
riod. 

We removed long-term trends, folded the APT pho¬ 
tometry over the orbital period (Porb = 19.6027 days), 
and binned the data (0.025 in phase). The resultant 
Johnson B and V light curves are presented in Figure 
O The quasi-sinusoidal trend observed in the averaged 
light curves suggests the possibility of ellipsoidal varia¬ 
tions due to distortions of the primary star partially fill¬ 
ing its Roche lobe potential. With a Roche lobe rad ius 
of 16.5 Rq, we obtain = 0.61 (|Eggletonlll983ll . 

We used the light-curve -fitting software package 
Eclipsing Light Gurve lELG: lOrosz fc Hauschildtll200(lll 
to model the ellipsoidal variations using our orbital 
parameters with no free param eters ( gravit y darken¬ 
ing assumed to be /3 = 0.08; iLucvl Il967l see Fig¬ 
ure 0. The characteristics of the ellipsoidal varia¬ 
tions with this model as compared to the light curve 
of cr Gem indicate that the long-term signature likely 
is indeed due to ellipsoidal variations, in contrast to 
previous suggestions that the periodicity at Porb/2 is 
due to active longitudes aligned with the orbit (e.g., 
Henry et al.lll99^ lJetsuiri996t iBerdvugina fc TuominenI 

1998HKaiatkari et al.ll20l4lKc)vari et al.ll2014[ l. We note 

that rotation pe riods derived from th e analysis of the 
light curve le.g.. IKaiatkari et al.ll2014ll suggest the star 
is rotating slightly faster than the orbital period, further 
supporting our identification of ellipsoidal variations in 
cr Gem. It should be noted that removing the effect of 
ellipsoidal variations from the light curve does not elim¬ 
inate all starspot signatures (See Appendix B). 

The ELG model fit of ellipsoidal variations can be im¬ 
proved to better match our data. We modeled the sys¬ 
tem again with no free parameters exc ept for the gravit y 
da rkening coefficient. 0 for T^ff oc (Ivon Zeipellll924ll . 
as lEsoinosa Lara fc Rieutordl (j2012l l recently suggested 
P ^ 0.21 for convective stars, sub stantially higher than 
the canonical f3 ^ 0.08 (jLucvlll9^ value assumed in our 
hxed-parameter fit. Although our average light curve is 
still contaminated by some residual spot modulation, we 
find that /3 = 0.02 ± 0.02 with error bars determined by 
bootstrapping over observing seasons of the 27 years of 
observation in the APT light curve. This value strongly 
rules out /3 > 0.1 for this system (see Eigure0. 

6 . GONCLUSIONS 

In this work, we have made the hrst visual detections 
of the secondary star of a Gem using interferometric and 
spectroscopic observations. We establish the first visual 
orbit by combining the interferometric detections with 
radial velocity data. The determination of orbital pa- 














































A V A B 


Detecting the Companion and Ellipsoidal Variations of cr Geminorum 


7 



Phase 

Fig. 5. — Differential folded and binned light curves of cr Gem 
for B and V magnitudes plotted with error bars from the binning. 
Each data point is an average of data points spanning 0.025 in 
phase from the complete light curve folded on the orbital period. 
The quasi-sinusoidal signature of the averaged light curve is due 
to ellipsoidal variations caused by the primary star partially filling 
its Roche lobe potential. The lines represents the ELC models 
for ellipsoidal variations with the gravity darkening coefficient /3 = 
0.02, 0.08, and 0.25, where /3 = 0.02 ± 0.02 is the best fit to the 
binned and averaged light curves. 

rameters has allowed for mo del-independent mass deter¬ 
minations (Ml = 1.28±0.07 Mq, M 2 = 0.73±0.03 Mq). 

Folded and binned photometric data have shown ev¬ 
idence of ellipsoidal variations, gravitational distortions 
of the primary star caused by the close companion. The 
light curve is comparable to light curve models created 
only from stellar and orbital parameters (assuming no 
starspots). Although the ellipsoidal variations are only 
a small effect, the primary star of a Gem is not spheri¬ 
cal, partially filling its Roche lobe potential and having 
a surface temperature gradient. Our establishment of el¬ 
lipsoidal variations offers a compelling alternative expla¬ 
nation to the previously purported detections of active 
longitudes, starspots on either side of the primary star 
(iHenrv et al.lll995HJetsulll996HBerdvugina fc TuoininenI 
119981 : iKaiatkari et al.ll20l4 iKovari et al.ll2014ll . 


Our new orbital elements along with the folded light 
curve also allow for measurements of gravity darken¬ 
ing. We find that [3 = 0.02 ± 0.02, a valu e of graviW 
darkening lower than suggested by th eory (|Lucvl HOCT 
lEspinosa Lara fc R.ieutordll201lL 120121 1. 

In this paper, we have demonstrated that precision in¬ 
terferometry at CHARA is now capable of detecting the 
faint main-sequence companions of bright RS CVn pri¬ 
mary stars. We are currently processing new data for 
other close, bright RS CVn systems and will be publish¬ 
ing these results in a series of follow-up papers. 
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APPENDIX 

A. CFA RADIAL VELOGITIES 

The CfA radial velocity data were obtained using two identical CfA Digital Speedometers (|Lathamlll992f) on two 
different telescopes: the 1.5-m Wyeth Reflector at the Oak Ridge Observatory located in the town of Harvard, Mas¬ 
sachusetts, and the 1.5-m Tillinghast Reflector at the Fred Lawrence Whipple Observatory (FLWO) on Mount Hopkins, 
Arizona. In the first set of radial velocity measurements, altogether 78 observations were obtained; 32 with the Till¬ 
ingha st Reflector and 46 with the Wyeth Reflector. Those radial velocities have all been published (iMassarotti et al.1 
EqoI), so the details of the procedures and reductions will not be repeated here, except to provide some overall char¬ 
acteristics of the data. Forty of the observations were obtained over a span of less than an hour using the Wyeth 
Reflector. They have been averaged to give one data point for the analysis in this paper. However, those forty obser¬ 
vations provide an opportunity to evaluate the precision of the individual velocity measurements from the CfA Digital 
Speedometers for an RS CVn primary with line broadening corresponding to a rotational velocity of 25 kms“^. The 
standard deviation of a single vel ocity from the average of all 40 is 0.40 kms“^. This compares favorably with the 
value of 0.45 kms“^reported bv IMassarotti et al.l (l2008f) for the RMS velocity residuals from their orbital solution 
using all 78 velocities. 

To avoid confusion, the velocities determined with the CfA Digital Speedometers have always been published on the 
native system of the instruments. For the analysis in this paper we added 0.14 kms“^to the published velocities, to 


^ Available at http://www.jmmc.fr/searchcal 


^ Available at http://cdsweb.u-strasbg.fr/ 
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put them on an absolute system defined by the lAU Radial-Velocity Standard Stars (iStefanik et al.lfloM . note that 
the sign of the correction was given as minus by mistake in that paper). 

A new set of velocities for both compon ents of the a Gem system was obtained with the Tillinghast Reflector 
Echelle Spectrograph (TRES; iFuresd[200811 at FLWO during the period 2012 October 1 to 2015 January 9. TRES 
is a modern fiber-fed CCD echelle spectrograph with resolution of 6.7 kms“^, very similar to that of the CfA Digital 
Speedometers. However, the free spectral range of the echelle order centered at 519 nm near the Mg b features is 
10 nm, compared to 4.5 nm for the CfA Digital Speedometers. For both instruments this is the wavelength window 
used for the determination of absolute velocities and for TODCOR analyses. Furthermore the signal-to-noise ratio 
per resolution element (SNRe) achievable with the TRES CCD detector is much higher than was possible with the 
intensified photon-counting detectors used in the CfA Digital Speedometers. The SNRe values for the old CfA spectra 
ranged from 30 to 100, while the typical value for the new TRES observations is 500. 

Fourteen strong TRES spectra were obtained of a Gem over a period 30 nights in 2014 December and 2015 January 
with the goal of detecting the lines of the secondary and deriving a double-lined spectroscopic orbit for the first time. 
Fortunately, two earlier observations from 2012 October were available in the TRES archive, which provided a two- 
year baseline for determ ining a more accurate orbital period. All 16 TRES spectra were analyzed using TODCOR 
(jZucker fc Mazehlll994D as implemented at CfA by G. Torres, and using the CfA library of synthetic spectra to choose 
the optimum templates. Only one of the observations failed to give a reliable velocity for the secondary, due to close 
blending of the lines from the two stars. 

The light of the primary dominates the composite spectrum of cr Gem, so it was straightforward to choose the 
template that gave the highest value for the average peak of the one-dimensional correlations; the parameters for the 
best template were Teff 4 =4500K,log5i = 2.5 (cgs),usinii = 25 kms“^, and solar metallicity. Because the secondary 
is so much fainter, its template parameters were only weakly constrained by the two-dimensional correlations. We 
therefore adopted the mass derived in this paper to guide the choice of secondary template parameters, Teff _2 = 4250 K, 
log 52 = 4.5 (cgs), usin *2 = 2 kms“^, and solar metallicity; the rotational velocity was selected under the assumption 
of tidal synchronization of the secondary spin with the orbital period. We also tried TODCOR solutions using the 
nearest neighbors from our library for the secondary template, and found that the results were not sensitive to the 
secondary template parameters. 

TODCOR has a mode in which the light ratio between the secondary and the primary in the observed wavelength 
window can be treated as a free parameter and thus can be determined from the spectra. We selected 10 observations 
where the velocities of the primary and secondary were well separated and used these to determine a light ratio at 
519 nm of 0.0139 ± 0.0025 (standard deviation of the values from the mean). If the errors are well behaved in this 
analysis, the uncertainty of the mean light ratio from TODCOR could be formally better by about a factor of 3, or 
about 6%. For the final TODCOR analysis we fixed the light ratio to 0.0139 for all the observations. The TODCOR 
velocities and times of observation are reported in Table|2]on the same zero point as the native CfA Digital Speedometer 
system. For the analysis in this paper, 0.14 kms“^was added to the velocities in Table [2] to put them on the lAU 
system. 

The orbital parameters are reported in Table [7] for two different double-lined spectroscopic orbits based on just 
the TODCOR velocities from the 16 TRES spectra. For the first solution the eccentricity was allowed to be a free 
parameter. Since the derived eccentricity was not significant, e = 0.0018 ± 0.0029, we also derived a solution for a 
circular orb it. Note that the mass r atio is well constrained; q = 0.582 ± 0.016 (2.7%). Table [7] also contains the orbit 
reported bv iMassarotti et ‘ahi (|2008f) and the orbit fit with only the data from the AST/TSU data set. Like the results 
from the TRES spectra, the eccentricity from the AST radial velocities, e = 0.002 ±0.003, is consistent with a circular 
orbit. Note that these are the entire data sets collected at these facilities and contain radial velocity measurements 
that were obtained when starspots were present on the stellar surface. 

Our double-lined orbital solution is among the most extreme that we have derived using the Mg b region in terms 
of the light ratio. Indeed, attempts to get good solutions from the neighboring echelle orders on either side of the 
Mg b order were unsatisfactory. Therefore we decided it would be prudent to test the reliability of our TODCOR 
analysis by creating sets of simulated observations using our library of synthetic spectra. For the time of each observed 
spectrum we simulated that observation by shifting the two synthetic template spectra by the velocities from the 
orbital solution and coadding after scaling by the light ratio found by TODCOR from the real observations. We found 
that a TODCOR analysis of the simulated observations reproduced the mass ratio and light ratio of the real data well 
within the estimated errors, and that the mass ratio was not sensitive to the value chosen for the light ratio. Thus 
we have no evidence for a systematic error in our determination of a light ratio of 0.0139 at 519 nm. This is quite 
different from the light ratio reported for the interferometric observations of 0.0038 at iJ-band. 

B. DIFFERENCE LIGHT CURVES 

In order to better justify our conclusion that ellipsoi dal variations can expl ain previous claims of “active longitudes” 
on a Gem, we have re-plotted some photometry from iKaiatkari et al.l ([20141 1 in Figure [6] along with our prediction of 
the expected ellipsoidal variation component using the ELC software and system parameters from Table|6]using gravity 
darkening parameter /? = 0.02. In FigurelHl we include data from two epochs , one showing very littl e overall variability 
and one showing high variability. In the first epoch (“Segment 8, Set 45” of [Kaiatkari et ahi (|2014f) L the photometric 
data showed clearly a d ouble-peaked light cur ve when phased with the orbital period, previously interpreted as due 
to active longitudes Isee lKaiatkari et al.ll2ni4l l. Here, we now see by removing the expected ellipsoidal variation, the 






















TABLE 7 

cr Gem Orbital Parameters from Radial Velocity Curves" 


Parameter 

TRES/TODCOR 
e free 

TRES/TODCOR 
e fixed 

Massarotti et al. f2008') 
e free 

AST/TSU 

e free 

AST/TSU 
e fixed 

orbital period, Porb (da^ys) 

19.6059 ± 0.0020 

19.6065 ± 0.0018 

19.60437 ± 0.00053 

19.6041 ± 0.0002 

19.6041 ± 0.0002 

center of mass velocity, 7 (kms“^) 

43.554 ± 0.077 

43.553 ± 0.074 

43.043 ± 0.066 

43.25 ±0.07 

43.33 ± 0.06 

semi-amplitude, primary, Ki (kms“^) 

35.19 ±0.10 

35.192 ± 0.098 

34.776 ±0.100 

34.60 ±0.08 

34.60 ± 0.08 

semi-amplitude, secondary, K 2 (kms~^) 

60.5 ± 1.6 

60.5 ± 1.5 




eccentricity, e 

0.0018 ± 0.0029 

0 

0.0143 ±0.0026 

0.002 ± 0.003 

0 

time of periastron passage, T (HJD) 

2456899.0 ± 5.9 


2453507.96 ±0.71 

2456985.3 ±3.7 


longitude of periastron, u)\ (°) 

38. ± 109 


46 ± 13 

181 ± 68 


time of maximum velocity To (HJD) 


2456916.571 ±0.012 



2453563.81 ±0.02 

Ml sin^ i {Mq) 

1.126 ±0.068 

1.126 ±0.065 




M 2 sin® i (Mq ) 

0.655 ±0.022 

0.655 ±0.021 




ai sin i ( 10 ® km) 

9.487 ±0.027 

9.488 ±0.027 

9.374 ± 0.024 

9.33 ±0.02 

9.33 ±0.02 

a 2 sin i ( 10 ® km) 

16.31 ±0.43 

16.32 ± 0.42 




a sin 2 (Rq) 

37.09 ± 0.62 

37.09 ± 0.60 




mass ratio, q = M 2 /Mi 

0.582 ±0.016 

0.582 ±0.15 




mass function, f(M) = (M 2 sin i)®/(Mi +M 2 )^ 



0.0854 ± 0.00066 

0.0841 ± 0.0006 

0.0841 ± 0.0006 

RMS velocity residuals, cri (kms~^) 

0.29 

0.28 

0.45 

0.30 

0.30 

RMS velocity residuals, (T 2 (kms~^) 

4.8 

4.6 




light ratio at 519 nm 

0.0139 

0.0139 





Note. — “These radial velocity data sets are the complete sets and include data contaminated by starspots. 
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Fig. 6 .— Average-sub tracted, differential lig ht curves of cr Gem for Johnson V magnitudes plotted for JD 2449802.6965 — 244985 9.6453 
(left; Segment 8, Set 45 JKaiatkari et alll2014l') ) and JD 2449982.9990 — 2450032.0349 (right; Segment 9, Set 1 l|Kaiatkari et alll2014l'l ). The 
top panel contains a plot of the APT data sets (circles) and the model ellipsoidal variations created with ELC for the orbital parameters 
of cr Gem and best-fit gravitational darkening coefficient /3 = 0.02 (solid line). The bottom panel contains the residuals of the APT light 
curve with the ellipsoidal variation signature removed (circles). 


TABLE 8 

Johnson Photometric Ellipsoidal Variations Models of a Gem 


Phase 

U 

B 

V 

R 

I 

J 

H 

K 

0.000 

-0.018 

-0.016 

-0.014 

-0.013 

-0.012 

-0.011 

-0.011 

-0.010 

0.011 

-0.019 

-0.016 

-0.014 

-0.013 

-0.012 

-0.011 

-0.011 

-0.010 

0.022 

-0.019 

-0.016 

-0.014 

-0.013 

-0.012 

-0.011 

-0.010 

-0.010 

0.033 

-0.018 

-0.015 

-0.014 

-0.012 

-0.011 

-0.011 

-0.010 

-0.009 

0.044 

-0.017 

-0.015 

-0.013 

-0.012 

-0.011 

-0.010 

-0.009 

-0.009 


Note. — Table is published in its entirety in the electronic edition. Or¬ 
bital phase is based upon our orbital parameters listed in Table Notably, 
To ^ 2453583.98 (HJD), Porb ^ 19.6027 days. 


signature of two spots on opposite sides of the star (the basis for th e active longitudes cla ims) nearly completely 
disappears (see Figure [6]). The second epoch (“Segment 9, Set 1” of iKaiatkari et all (|2014lb is dominated by one 
spot and the ellipsoidal variations are not discernible. Nonetheless, future starspot modelers should account for the 
underlying ellipsoidal variations before performing detailed light curve analysis or surface brightness inversions. We 
have included our calculation of the expected ellipsoidal variations for Johnson UBVRIJHK in Table |5] to assist 
future workers-a full re-analysis of the star spot properties is beyond the scope of this paper. 
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